Numerical Contractor Renormalization Method for Quantum Spin Models

We demonstrate the utility of the numerical Contractor Renormalization (CORE) method for quantum spin systems by studying one and two dimensional model cases. Our approach consists of two steps: (i) building an effective Hamiltonian with longer ranged interactions up to a certain cut-off using the CORE algorithm and (ii) solving this new model numerically on finite clusters by exact diagonalization and performing finite-size extrapolations to obtain results in the thermodynamic limit. This approach, giving complementary information to analytical treatments of the CORE Hamiltonian, can be used as a semi-quantitative numerical method. For ladder type geometries, we explicitely check the accuracy of the effective models by increasing the range of the effective interactions until reaching convergence. Our results in the perturbative regime and also away from it are in good agreement with previously established results. In two dimensions we consider the plaquette lattice and the kagomé lattice as non-trivial test cases for the numerical CORE method. As it becomes more difficult to extend the range of the effective interactions in two dimensions, we propose diagnostic tools (such as the density matrix of the local building block) to ascertain the validity of the basis truncation. On the plaquette lattice we have an excellent description of the system in both the disordered and the ordered phases, thereby showing that the CORE method is able to resolve quantum phase transitions. On the kagomé lattice we find that the previously proposed twofold degenerate S = 1/2 basis can account for a large number of phenomena of the spin 1/2 kagomé system. For spin 3/2 however this basis does not seem to be sufficient anymore. In general we are able to simulate system sizes which correspond to an 8 × 8 lattice for the plaquette lattice or a 48-site kagomé lattice, which are beyond the possibilities of a standard exact diagonalization approach.